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ABSTRACT 

Atmospheric optical turbulence induces phase fluctuations in a propagating elec- 
tromagnetic wave. The resulting degradation in coherence limits the capability of any la- 
ser, target acquisition, or surveillance system. Past data collection methods for the 
parameterization of atmospheric turbulence profiles, in support of critical Theater Ballistic 
Missile Defense (TBM) systems, from ground level to 30 km, have depended on mete- 
orological balloon-thermosonde systems, probes carried on the U.S. Air Force Argus air- 
craft, as well as radar and optical measurements.. The balloon and aircraft systems 
measure temperature fluctuations and compute the temperature structure function, Ct 2 and 
the related index of refraction structure parameter, C n 2 . It has recently become critical to 
explain why turbulence profiles from daytime thermosonde data consistently show a two 
order of magnitude increase over that taken during the night, primarily between 12-20 
km. 

This thesis analyzed the TSI 3.8 pm platinum coated tungsten thermosonde probe 
used by the USAF Research Laboratory (AFRL) to quantify the magnitude of the solar 
heating and to investigate other heat transfer mechanisms in the probe. A model of the thin 
wire probe was developed to identify each of the contributions to the temperature error 
and its significance. Experimental measurements where collected to verify most aspects of 
the final model. 

We found that the sun induces a temperature rise in the TSI 3.8 pm fine wire 
probe, during the day, that can vary from near zero to 0. 175 K. It is strongly dependent on 
probe orientation with respect to the sun and on variations in the air flow over the probe. 

This then causes an apparent increase by two orders of magnitude in the daytime meas- 
urements of the optical turbulence parameters Ct 2 and C n 2 
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I. 



INTRODUCTION 



The Airborne Laser Program (ABL) must contend with several problems associ- 
ated with the propagation of electromagnetic radiation through the atmosphere. Specifi- 
cally, random fluctuations in the atmospheric index of refraction cause phase front 
distortions of the propagating plane wave [Ref. l:p. 6-6], Other effects such as thermal 
blooming, absorption and scattering by atmospheric particles, and stimulated Raman scat- 
tering are also important but go beyond the scope of this thesis [Ref. 2], The concept of 
mounting a weapons class laser on an aircraft has been around for over two decades and 
was proven possible by the Airborne Laser Laboratory in the 1970’s. The technologies to 
make it practical, however, have only recently been demonstrated by Lockheed Martin and 
scale laser testing done by the Air Force Research Laboratory in 1995 as part of the 
ABLE ACE experiment [Ref. 3:p. 6], 

The basic requirements of the ABL system are to detect, track, and destroy theater 
ballistic missiles (TBMs) while still in their boost phase and prior to the missile’s leaving 
the hostile country’s own airspace. The ABL system will use a passive infrared (IR) 
search and track system along with an active Multi-Beam Illuminator (MBEL) laser for 
accurate positioning and better dwell characteristics during high power Chemical Oxygen 
Iodine Laser (COIL) missile destruction [Ref. 4:p. 1]. Adaptive optics will be needed to 
achieve effective ranges beyond ~ 100km. 

The ABL program utilizes atmospheric turbulence data that have been collected 
by the Air Force Research Laboratory (AFRL), Naval Postgraduate School (NPS), and 
Argus aircraft for this and other Department of Defense atmospheric laser programs over 
the past decade. Characterization of atmospheric turbulence is essential in the design and 
use of atmospheric laser systems. A critical factor in the design is whether a perceived 
enhancement in the measured temperature structure parameter, Cj 2 , during the daytime as 
compared to night is real or is a consequence of solar heating of the measurement appa- 
ratus itself [Ref. 5:p. 30], Experiments done by AFRL, as well as subsequent ones done 
by NPS, have shown a significant enhancement in daytime Ct 2 values, most likely caused 
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solar heating of the thermosonde probes used to do the measurement, coupled with the 
complications of package rotation, pendulum motion, and probe vibration [Ref. 5:p. 55]. 

This thesis investigated the systematic effects caused by solar heating of the bal- 
loon thermosonde probes used by AFRL and NPS. Extensive probe response modeling 
and experimental measurements were done to identify and to quantify the individual con- 
tributions to temperature variation and to determine if the signal generated would be in- 
terpreted as a false turbulence signal within the electronic bandpass of the instrumentation. 
Package rotation, pendulum motion, and other dynamics where apparent in the experi- 
mental data. 

The ABL program has entered the Program Definition and Risk Reduction phase 
(PDRR) which calls for one complete weapon system, including a laser and its controlling 
system, to be integrated onboard a 747-400 platform and to be successfully demonstrated 
by 2002. Success of the system will be evaluated by attempting to shoot down a SCUD 
type missile [Ref. 6:p. 1-2]. Essential risk assessment and management issues underlying 
the feasibility of proceeding past the project’s scheduled acquisition program Authority to 
Proceed Milestone (ATP) in mid- 1998 relies partly on the resolution of the atmospheric 
turbulence problem. The underlying motivation of this investigation was to identify any 
systematic error sources that might affect the microthermal platinum coated tungsten 
probes that provide critical turbulence measurements needed for effective deployment of 
theater ballistic missile defense systems. 
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H. BACKGROUND 



A. OPTICAL TURBULENCE PARAMETERS 



Atmospheric turbulence degrades the wave front of a collimated laser beam Ran- 
dom fluctuations in the atmospheric index of refraction, n, caused by temperature and 
humidity inhomogeneities, induce centroid wander, beam broadening, and interference 
that all reduce the beam’s spatial coherence (Figure 2-1) [Ref. 7:p. 1], 

Tatarski and others have developed useful parameters for the statistical quantifi- 
cation of optical turbulence in the atmosphere by defining so-called structure functions. 
By describing the atmosphere as a three dimensional random field Tatarski considered the 
correlation function of some relevant quantity of interest, f, as [Ref. 7:p. 3] 

)>][/ ‘ (< 2 )-</■(>,))]) . ( i ) 

where the angled brackets denote the mean value, the superscript star is the complex con- 
jugate, and t denotes time. Because the atmosphere is not statistically homogenous or 

stationary, the mean values (/(*,)) are not well defined. To avoid this problem Kol- 
mogorov and Tatarski introduced the structure function 

D, (<„<,)={[/(<, )-/(/,)]’}. (2) 

which characterizes the intensity of the particular random fluctuations of f(t) [Ref. 7:p. 9- 
10], For our use then, a structure function tensor can be described by the mean square 
difference between two measurements, such as temperature, by 

D T (r„r,) = (imr,)-T(r 2 )f). (3) 
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Figure 2-1: Effect of atmospheric turbulence on a coherent plane wave. 



For an isotropic and homogeneous quantity the structure function must necessarily be in- 
dependent of the orientation angle and absolute position, and can depend only on r, the 
distance between the two observation points: 

r = |F— r 2 |. (4) 

Although the atmosphere is generally not isotropic or homogenous, if r is sufficiently 
small, then some local volume may be sufficiently isotropic and homogenous to permit 
credible measurements [Ref. 7:p. 31]. 

It is found that the statistical behavior of the atmosphere is well behaved over a 
limited range of r, l 0 <r < L 0 . This is called the inertial subrange. The inertial subrange 
is bounded on the lower end by / 0 , which is the small scale turbulence region associated 
with viscous eddy dissipation. On the upper end of the inertial subrange, Lo defines the 
outer scale length, which is the largest eddy size that can be considered isotropic through- 
out its volume [Ref. 8], 

Kolmogorov and Obukhov proposed that the structure functions in the atmosphere 
have a r 2/3 dependence [Ref. 7:p. 31]. Assuming isotropic and homogenous conditions, 
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Tatarski proposed the temperature structure function, Cj 2 by normalizing Eq. (3) by r 23 



[Ref. 7:p. 46] 



p T (f) 

L*T ~~ 2 2 



( 5 ) 



r 



Since n - 1 is proportional to density, which is a function of the pressure, P, in 
millibars and temperature, T, in K, the inertial range refractive index structure parameter, 



where z is the altitude in meters [Ref. 7:p. 79], The temperature structure parameter, Ct 2 , 
is measured by fine wire temperature sensors in a balloon or aircraft system as a function 
of altitude. C n 2 is then computed from Ct 2 using Eq. (6). 

The effects of turbulence on a laser beam propagating through the atmosphere ap- 
pear as distortion in the final beam pattern at the target. A common measure of beam 
quality is the Strehl ratio, which is the relative intensity of the beam’s point spread func- 
tion measured at the mean centroid location relative to that of an ideal, diffraction limited 
optical system [Ref 9], The Strehl ratio depends on the spatial coherence of the laser ra- 
diation after passing through the atmosphere, which in turn depends on the atmospheric 
turbulence and its effects on the modulation transfer function [Ref. 18], The atmospheric 
modulation transfer function (MTF), characterized by the transverse coherent path length 
t 0 , is a measure of electric field transverse autocorrelation length, [Ref. l:p. 6-29]. Fried 
showed that, for plane waves. 



C n 2 , is related to the temperature structure parameter by 




( 6 ) 



L 




(V 



0 



where k is the optical wave number (2 tc/X.) and L is the optical path length [Ref. 10:p. 
397], Values of r 0 range from centimeters under very turbulent conditions to meters 




Figure 2-2: Physical representation of the coherence length, r 0> and its reduction due to 
atmospheric turbulence. 

when turbulence is very low. Since r 0 describes the spatial coherence of the beam, the 
larger its value the better (Figure 2-2). The diffraction limit of the system, however, de- 
pends on the geometry of the optical system as well as r<j. Without turbulence the diffrac- 
tion-limited resolution in radians for a circular aperture is [Ref. 1 1 :p. 393] 

1 . 222 . 

. ( 8 ) 

where X is the wavelength and D is the aperture diameter. When D>r 0 , the effective tur- 
bulence-limited angular resolution just becomes 

A (9) 

r o 

The isoplanatic angle, 0 O , is another useful statistical parameter for the characteri- 
zation of atmospheric optical turbulence. It is related to r 0 , and represents a measure of 
the angular field of view of an optical system. It characterizes a cone within the atmos- 



6 



phere with similar turbulence characteristics. Fried [Ref. 12] showed that 0 O depends on 
the integration of C n 2 along an optical path and is given by 

r L I- 3 ' 5 



*0 = 



2.905 k 2 jc 2 n (z)z $,3 dz 



( 10 ) 



with L, the optical path length in meters. 

As has been previously stated, current balloon and aircraft probe thermosonde 
systems can measure Ct 2 and then compute C„ 2 as a function of altitude and time. Using 
Eqs. (8) and (9) this data provides the basis for computing r 0 and 0 O for a particular appli- 
cation in the design and development of an adaptive optics system. 



B. THERMOSONDE 



The AFRL thermosonde design used during the last decade consists of a VIZ 8000 
electronics package with two fine wire TSI 3.8 pm platinum coated tungsten probes ex- 
posed to the ambient atmosphere and separated by 1 m (Figure 2-3). This system meas- 
ures the wires’ resistance fluctuations caused by changes in temperature. The balloon 
carrying the package is capable of ascending to an altitude of approximately 30 km. 

Changes in the probe resistance, Rw, give the temperature, Tw, through the rela- 
tion 



RfV - ^0 + [\ + a (T iV ~ T 0 )] 



( 11 ) 



[Ref. 13:p. 11], where Ro and To, are values at a specific reference resistance and tem- 
perature respectively and a is a constant for the particular wire composition and geome- 
try, platinum coated tungsten in this instance. In practice, the two probes form the 
opposing elements of an AC Wheatstone bridge. The difference in resistance from the 
Wheatstone bridge is amplified and converted to a temperature change by [Ref. 14:p. 15] 



AT = 



oR o 



( 12 ) 
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The probe data channels are sampled every 10 s along with ambient pressure, mean tem- 
perature, and humidity. The data are radioed via the modified VIZ 8000 radiosonde to the 
ground station that continuously displays and stores the data on a PC. Optimum ther- 
mosonde performance requires a compromise in the probe current. A value of 0. 1 ma for 
a 20 Cl probe provides sufficient signal without inducing temperature differences above 
0.002 K between the probe and the atmosphere. The ohmic heating temperature differ- 
ence, ts.T w , is 



A T w 



q 



(13) 



where q w is the electrical power input, As is the surface area of the wire, and h is the coef- 
ficient of heat transfer [Ref. 14:p. 17], 

The response of the probe to a change in temperature must be faster than the ac- 
tual fluctuation to be measured if a meaningful measurement is to result. Failure of the 
probe response to follow the actual atmosphere will cause an underestimate of the struc- 
ture function. The response time of our probe is given as [Ref. 14 :p. 20] 



T - 



r nD 2 L? 
4 > 

h{7dDL) 



pDc 
~ 4 ~h 



(14) 



with: 

p = density of probe wire (kg/m 3 ), 
c = heat capacity of wire (J/kg K, 
h = convective heat transfer coefficient (W/m 2 K), 
D = wire diameter (m). 
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Figure 2-3: Construction layout of platinum coated tungsten AFRL probe and support. 

C. CURRENT APPLICATIONS 

The apparent increase in turbulence during the day compared to night has signifi- 
cant implications on the ABL program (Appendix A). Current performance specifications 
for the ABL’s COIL laser require the destruction of a ballistic missile that is still in it’s 
boost phase out to a range of several hundred kilometers [Ref. 15:p. 1] The power re- 
quirements on target for these ranges have only recently become possible with the use of 
adaptive optical systems that rely on active compensation of wave front distortion of the 
propagating laser beam. Still, only under favorable atmospheric conditions can ranges 
such as these be attained. If existing daytime balloon measurements of Cj 2 are valid they 
imply that the useful range of the laser system is severely limited in the daytime. Other new 
programs which rely on laser propagation through the atmosphere to great distances 
would also have to be substantially re-thought. 



9 



D. PREVIOUS RESULTS 



Numerous experiments have been performed by AFRL and NPS to characterize 
the microthermal probes and to determine Ct 2 and C n 2 during day and night conditions. 
Brown has flown many balloon launches since July of 1981 and performed extensive 
analysis to show that the probe supports do not contribute to observed diurnal differences 
[Ref. 5, 14, and 16], Walters, Weitekamp, Roper, and others through multiple balloon 
launch experiments (notably from 31 May to 3 June 1990 at AMOS observatory on Mt. 
Haleakala, Hawaii and 4 -9 December 1990 at Holloman Air Force Base, New Mexico) 
showed that solar contamination of data was a real effect. Figure 2-4 is a representative 
example of night time data collected over Maui on 2 June 1990 showing that at least some 
of the raw C n 2 data reached the lower noise limit (solid line) of the sonde electronics. Fig- 
ure 2-5, shows a two order of magnitude shift in C n 2 data for a next day launch from the 
same location. None of the data approached the noise limit of the electronics. One would 
expect that under normal conditions, there should be some parcel of air within the atmos- 
phere that is not turbulent. Yet the measurements show that that assumption is valid at 
night but not during the day. Hence some doubt has been cast on the general validity of 
the daytime measurements. 

Additional work by Walters, Weitekamp, and Roper have also shown that balloon 
wake effects are not a source of contamination [Ref. 17 and 18]. Additional attempts at 
methods of sonde stabilization have also been investigated by them to eliminate solar ef- 
fects, but to date no satisfactory method has been found. 
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Figure 2-4: Raw night time AF'RL probe C„ 2 data collected on Maui, 2 June 1990 show- 
ing the lower electronic noise limit (solid line) and corresponding data points near this. 
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Figure 2-5: Raw day time AJFRL probe C„ 2 data collected on Maui, 2 June 1990 showing 
the lower electronic noise limit (solid line) and corresponding data points with solar in- 
duced shift of two orders of magnitude. 
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m. MODELING THERMOSONDE PROBE HEAT TRANSFER EFFECTS 



A. INTRODUCTION 

In this Chapter, the heat balance equation for the TSI probe wire is described in 
detail. Special attention is paid to convective cooling of the wire, the effects of changing 
solar illumination areas, and effective free stream velocity fluctuations. MATLAB pro- 
gramming was used to model the probe under these changing conditions as a function of 
altitude. Once the model was validated by experimental measurement as described in 
Chapter IV, simulations of different probes’ responses to material and geometry changes 
was investigated as shown in Chapter V. 

B. HEAT BALANCE EQUATION 

The equilibrium temperature of a fine wire probe is given by the heat balance 
equation 



where: 

Qconv = convective transfer of power from the probe to the atmosphere (W), 

Q solar = solar power absorbed directly by the probe (W), 

Qrefi = reflected solar power from the Earth absorbed by the probe (W), 

Qpianck = Planck radiation power transfer between the probe and the environment (W), 
Q ohmic = ohmic heating caused by current in the probe wire (W). 

In terms of power (W) the equilibrium condition is 



Qconv Q solar Qrefi Qpianck Ql 



Ohmic 



= 0 , 



( 15 ) 




( 16 ) 
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The individual terms in Eq. (16) will be defined in the section in which the heat 
transfer mechanism is introduced. The thin wire probe areas Ai -A$ depend on the heat 
transfer mechanism, package pendulum motion and rotation, probe wire vibration, and 
solar zenith angle and are also defined in their appropriate sections. 

The time for wire to reach thermodynamic equilibrium is short compared to the 
heat transfer time via convection, so equilibrium was assumed [Ref. 17:p. 32], Thermal 
conduction through the probe wire to the supports was also ignored as shown in [Ref. 

5:p. 45] for the simulation, but was experimentally measured as described in Chapter IV. 

During flight, the effective free stream air velocity over the probe wire is also vari- 
able. Variations in the wind velocity alter the heat transfer from a probe that has a different 
temperature from the air. These fluctuations were modeled both as a steady state charac- 
teristic air velocity (V= 5 m/s) and as a random fluctuation distribution about a character- 
istic free stream velocity (Figure 3-1). 




-4 -2 0 2 4 6 8 10 12 14 



Air Velocity (m/s) 

Figure 3-1 : Simulation of air velocity fluctuation over probe 
wire for a fluctuation of 2.5 m/s about a velocity of 5 m/s. 
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c. 



CONVECTIVE COOLING 



1. Convection with Steady State Free Stream Velocity 

The convective transfer of power from the probe wire to the atmosphere is 



Qconv=hAiAT, 



( 17 ) 



where: 



h = convective heat transfer coefficient (W/m 2 K), 

A x = 7zDL = surface area of the probe subject to convection (m 2 ), 
AT=T probe -T (K). 



Heat transfer from a heated probe wire to the atmosphere consists of free convection from 
buoyancy effects at low air velocity, forced convection at higher air velocities (Reynolds 
numbers. Re >0.5), and mixed free and forced convection between both these regions, 
about Re=0.5 [Ref. 20:p. 323], Figure 3-2, developed by Collis and Williams, illustrates 
the complicated behavior of this low Re region. 

The Reynolds number is the ratio of inertial to viscous forces on a parcel of fluid. 
In this application. 



where D is the probe diameter, v, is the kinematic viscosity (which is the ratio of dy- 
namic viscosity, p, and density, p) and V is the characteristic velocity [Ref. 21:p.54], The 
NPS-AFRL probe, as modeled, would experience Reynolds number as low as 0.00024 but 
no more than about 4. The Nusselt number 



is another useful dimensionless number describing the heat transfer coefficient, where k t is 
the thermal'conductivity of the fluid, in this case air. 



DV 



(18) 



Re = 



v 




(19) 
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Figure 3-2: Interaction of free and forced convection from 
[Ref. 27], 
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At low Reynolds numbers, where buoyancy becomes important, the dimensionless 
Grassof number, given by 



Gr = 



P 2 sPziT probe - T)D 3 

2 



( 20 ) 



relates the ratio of buoyant to viscous forces near the wire where g is the local gravita- 
tional acceleration and p 2 is the coefficient of expansion of the fluid. The Grassof number 
is important when Re<Gr 1/3 and buoyant effects become a significant contributor to the 
probe’s ability to dissipate heat [Ref. 22:p. 131], Grassof numbers modeled for our probe 
varied from near zero to 3. 25x1 O' 8 . The effect of the Grassof number on the buoyancy ef- 
fect for low Reynolds number is also shown in Figure 3-2. At altitudes above 30 km mo- 
lecular diffusion becomes dominant over convective effects for the heat transfer 
mechanism of the AFRL probe wire [Ref. 14:p. 19]. 



2. Determination of h 

The convective heat transfer coefficient h is the heat flux leaving the wire surface 
per unit temperature difference between the wire and the free stream. The effect of h on 
the probe is significant and is a complex function of geometry, Reynolds, Grassof, and 
Prandtl numbers as well as the overheat ratio, (T OTO be-T)/T [Ref. 23: p.352], The dimen- 
sionless Prandtl number, Pr, is the ratio of the momentum diffusivity to the heat diffltsivity 
and is given by 

Pr= ir (21) 

The conductivity k t is defined in Table 3-1, and Cp is the specific heat of the atmosphere. 

The convective heat transfer coefficient can be found using a variety of methods. 
Most relate the Nusselt , Reynolds, and Prandl numbers and are valid over certain Rey- 
nolds and Nusselt number regimes, probe geometries and fluid/gas compositions. No sin- 
gle credible model for our probe’s entire operating environment was found. Four different 
models were compared, and the relative merits of each were contrasted to our probe’s 
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actual performance. Consideration was also given to the physical processes that were in- 
cluded in the model. The four methods considered are summarized below. 

a. h, Using Kothandaraman and Subramanyan’s Method 

Kothandaraman and Subramanyan in [Ref. 24: p. 91] proposed the follow- 
ing relation for the Nusselt number 

Nu D =-^- = CRe n , (22) 

where k t is the conductivity of the fluid (air) determined empirically in Table 3-1: 



Altitude (km) 


0-12 


12-20 


20-32 


k t (W/mK) 


Alt- 50.57 
- 1983.47 


0.0195 


Alt + 2065 

+ 11636 



Table 3- 1 : Conductivity k t of atmosphere variation with altitude in kilometers, adapted 
from the Geophysics Handbook [Ref. 25:p. 14-6], 



C and n are constants as given in [Ref. 24:p. 91] for the Reynolds number range 
0.4<Re<4.0 and have values of C=0.891 and n=0.330. Values of Re less than 0.4 require 
further extrapolation, as done in [Ref. 17]. The altitude, Alt, here and throughout, is 
given in kilometers. A characteristic velocity V of 5 m/s, corresponding to the approxi- 
mate ascent rate of the balloon, and characteristic length D of 4.7 pm (measured fine wire 
probe diameter; refer to Chapter V) were used. A more elaborate treatment for modeling 
the effective velocity over the probe wire was done in Ref. 5:p. 7-9 and Ref. 19:p. 22, but 
for our methods this characteristic velocity will be sufficient. The density assumed the va- 
lidity of the ideal gas law 



P = 



pRT 
M ’ 



(23) 



and modeled the pressure variation with altitude by 

P = p 0 e". 



(24) 
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R, above, is the universal gas constant (8.31 J/Mole-K) and M is the mean molecular 
weight of the atmosphere (28.9644 kg/kMol) [Ref. 25:p. 14-5], Atmospheric tempera- 
ture variation with altitude, T, came from the USAF Geophysics Handbook and is shown 
in Table 3-2. The scale height, zo, and scale pressure, p 0 , are 8.1 km and 1013 mBar 
(101.3 kPa) respectively [Ref. 25:p. 14-5], 



Altitude (km) 


0-11 


11-20 


20-32 


T (K) 


Alt- 44.62 


217 


Alt + 142.75 


-0.15493 


0.75 



Table 3-2: Atmospheric mean temperature variation with altitude adapted from the Geo- 
physics Handbook [Ref. 25: p. 14-3], 



Dynamic viscosity was determined using the relation 

3 



= M1 

M T + S ’ 



(25) 



[Ref. 25 :p. 14-5], Here pi is a constant (1.45x10^ kg/s-m-K 1/2 ) and S is Sutherland’s 
constant (1 10.4 K). 

With all quantities known and solving for h we find 



h = 



CRe” k 

D 



(26) 



For this probe the diameter, D, was 4.7 pm and the length, L, was 4. 15 mm. 



b. h, Using Kramer’s Method 

Kramer’s law for circular cylinders of infinite length is given by Hinze in 
[Ref. 26:p. 201] as 

Nu = 0.42 Pr 0 ' 20 + 0.57 Pr° 33 Re 0 50 . (27) 

Kramer’s law is valid for Reynolds numbers in the range 0.1<Re<10,000 [Ref. 19:p. 63], 
So now using Eq. (23), h can be found. 
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c. h, Using Zijnen’s Method 

The heat transfer proposed by Zijnen does not have any stated valid Rey- 



nolds number range for circular cylinders and is given as [Ref 19:p. 64] 



Nu = 0.3 8 Pr 0 2 + (0.5 6 Re 0 5 + 0.00 1 Re) Pr 0 333 . 



(28) 



Once again Eq. (23) is used to solve for h. 



d. h, Using Collis and William’s Method 

Collis and Williams developed the following relation for the Nusselt num- 



ber [Ref. 27:p. 357-384] 




(29) 



For air and in the Reynolds number range 0.02<Re<44 the constants are as follows: 
n=0.45, A=0.24, and B=0.56 [Ref. 22:p. 130]. The subscript f denotes measurements 
were taken at the film temperature which is the average temperature between the wire and 
the atmosphere. Eq. (23) then gives the relation for h in this method. 

Bailie further modified this relation to account for the mixed convection effects in 
the region around Re<0. 15, as shown in Figure 3-2, by introducing an effective Reynolds 
number, Re ef f, given by 



Bailie also investigated the error associated with assuming an infinite cylindrical probe in 
the above approximations and developed Figure 3-3, which shows the effect of aspect ra- 
tio (21/D) on the Nusselt number [Ref. 22:p. 132], The AFRL probe used has an aspect 
ratio of over 2000 and therefore behaves as an infinite wire to high degree. 



Re] ff = Re} + (0.9 Gr 0418 ) 2 . 



(30) 
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Figure 3-3: Effect of wire aspect ratio on Nusselt 
number. From [Ref. 22:p. 132], 

e. Comparison of h Values 

Although no method covered the extremely low Reynolds numbers en- 
countered in this work, the Collis and Williams method best represented the physical ef- 
fects encountered and the behavior of our probe that is as most represented by the low 
Reynolds number regions shown in Figure 3-2. All four methods were in good agreement 
at Reynolds numbers greater than 1. Final values for the coefficient of heat transfer for air 
velocities of 0 and 5 m/s as derived from the Collis and Williams technique are shown in 
Figure 3-4. 
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The kink in the curves in Figures 3-4 to 3-9 around 12 km comes from the rate of 
change of temperature versus altitude of this part of the stratosphere. 

Figure 3-5 summarizes the coefficient of heat transfer values for the four methods 
as a function of altitude for an air velocity of 5 m/s using a MATLAB program (Appendix 
B). This program was also used to produce the plots of the. Re, Nu, and Gr numbers, at- 
mospheric temperature T and kinematic viscosity v versus altitude, as shown as Figures 3- 
6 to 3-1 1 for the steady state free stream velocity of 5 m/s and as a function of altitude. 

Although not a factor with balloon mounted probes, those mounted on aircraft, 
such as Argus, also need to consider the heat loss effects at different Mach numbers. 
Spangeuberg [Ref. 22:p. 133] presents the reported heat loss at various Mach numbers 
and is reproduced here as Figure 3-12. 




Figure 3-4: Convective heat transfer coefficient of 4.7 pm probe 
using Collis and Williams method versus altitude with constant 
free stream velocity of 0 and 5 m/s. 



22 



Altitude (km) Altitude (km) 




Figure 3-5: Comparison of h, for the 4.7 pm probe, using 
the four methods investigated at an air velocity of 5 m/s. 




Figure 3-6. Reynolds number versus altitude, for the 4.7 pm probe, 
for a constant free stream velocity of 5 m/s. 
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Figure 3-7: Reynolds number versus altitude, for the 4.7 (im 
probe, for an air velocity of 0.01 m/s. 




Figure 3-8: Nusselt number versus altitude, for the 4.7 pm probe, 
for a constant free stream velocity of 5 m/s. 



24 




30 



25 




Figure 3-9: Grassof number, showing ratio of buoyancy to vis- 
cous forces, versus altitude for constant free stream velocity of 
5 m/s and corresponding to a differential probe temperature of 
0.035 K for the 4.7 pm probe. 




Figure 3-10: Atmospheric temperature versus altitude. 
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Figure 3-11: Kinematic viscosity of the atmosphere versus altitude. 




Figure 3-12: Heat loss at increasing Mach number, from Davies and 
Fischer in [Ref. 22 :p. 133], 
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3. Convective Cooling Power 

Typical values for the convective cooling of the wire at sea level are on the order 
of 10. 63x10^ W, but they are sensitive to the free stream velocity. Shown as a function of 
altitude in Figure 3-13 are values of Q con v for a free stream velocity of 5m/s. 




Figure 3-13: Probe convection, Q con v versus altitude (V = 5 m/s) 

4. Convection with Normal Free Stream Velocity Fluctuation 

Real world conditions encountered by the probe involve dynamic fluctuations of 
the free stream velocity component over the probe wire. These fluctuations arise from 
package pendulum motion, rotation, possible balloon wake turbulence, actual atmospheric 
turbulence, and to a lesser extent probe vibration. 

Velocity fluctuations ultimately affect the coefficient of heat transfer and Reynolds 
number which both depend on air velocity over the probe wire. The effect of a random 
normal mean free stream velocity fluctuation of 2 m/s about 5 m/s is shown in Figures 3- 
14 and 3-15 for both these quantities. The convection, Q CO nv, term will vary similarly with 
altitude and is not plotted here. 
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Figure 3-14: Modeled convective heat transfer coefficient with 
velocity variance of 2 m/s about 5 m/s for the 4.7 txm probe. 
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Figure 3-15: Modeled Reynolds number with normal mean ve- 
locity variance of 2 m/s about 5 m/s for the 4.7 pm probe. 
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D. DIRECT SOLAR ABSORPTION OF ENERGY 



1. Constant Probe Cross Sectional Area Presented to Solar Radiation 

The power absorbed directly from the sun by the probe wire is given by 

Qsolar=aA 2 P s , (31) 

in the heat balance equation with: 

a= solar absorptivity of the probe, 

A 2 = cross sectional area of probe seen by the sun (m 2 ), 

P s = solar power flux (W/m 2 ). 

The solar absorptivity of the platinum coated tungsten wire probe is 0.435 from tabulated 
values of the normal incident reflectivity (R=0.565 @ 3.0eV) for platinum wire, which is 
related to the absorptivity by a=T-R [Ref 28:p. 10-282], The maximum illuminated area 
of the probe as seen by the sun to good approximation is the cross sectional area DL. The 
projected area, A 2 , is a function of the complex probe pendulum motion and rotation, 
probe vibration, and solar zenith angle. The change in area subjected to solar radiation, as 
a function of time, was modeled as one half the magnitude of the sum of two sinusoidal 
variations of the total visible surface area. Periods for the two sine functions correspond 
to the observed frequency of pendulum motion and rotation. Theta dependence (sun ze- 
nith angle) was not modeled explicitly in this thesis but was done previously [Ref. 5:p. 10]. 
The solar flux, P s> varies with altitude (meters) as given in [Ref. 25] 

P s =(Alt) 05 * 4 +890, (32) 

to our maximum altitude of interest, 30 km. 

Typical values of the power absorbed directly from the sun by the probe at sea 
level are of the order 7. 28x10“* W and are shown as a function of altitude in Figure 3-16. 
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Figure 3-16: Direct solar absorption, Q so ur, of power with constant 
exposed cross sectional area versus altitude for the 4.7 pm probe. 

2. Variation of Probe Cross Sectional Area Presented to Solar Radiation 

The probe cross sectional area exposed to solar illumination is a function of the 
roughly sinusoidal geometry variations caused by pendulum motion and rotation with re- 
spect to the sun. Secondary influences of probe wire vibration are discussed below. Us- 
ing MATLAB programs developed by Walters (Appendix C), digital filters and Fourier 
spectral techniques were used to analyze single probe thermistor data taken at Anderson 
Peak in September 1992. This launch was during the day and showed the characteristic 
C„ 2 enhancement obtained from higher order structure functions (Figure 3-17). Figure 3- 
1 8 shows the resultant power spectrum plot with two distinct periodicities of about 5 and 
20 seconds. Visual observation of these balloon launches and simple calculation of the 
pendulum motion confirm that these correspond to the periods of the rotational and pen- 
dulum motion of the sonde package during its ascent. Hence these were used in our 
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Magnitude 



model. The influence of the sinusoidal variation of the solar cross sectional area A 2 on 
solar absorption is shown in Figure 3-19. 




Figure 3-17: Anderson Peak C n 2 during the day obtained 
from a single flake thermistor. 




Figure 3-18: Power spectral plot of Anderson Peak data 
- from Figure 3-17. 
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Figure 3-19: Effect of Sinusoidal variation of probe solar illumina- 
tion area, A 2 on solar heat flux with air velocity of 5 m/s for the 4.7 
|im probe. 

3. Probe Vibration 

Less well understood is the effect on solar heating caused by the probe wire vi- 
brating in and out of solar view. Experimental work done to measure this phenomenon is 
discussed further in Chapter IV. 




E. REFLECTED SOLAR ABSORPTION 



The reflected solar power from the Earth onto the probe wire is given by 

Q re fi=aA4PsY, (33) 

where: 

a= solar absorptivity of the probe, 

A 4 = effective area of the probe that absorbes reflected sunlight (m 2 ), 

P s = solar energy flux (W/m 2 ), 

y= percent of solar energy reflected by the Earth (terrestrial albedo) 
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The solar absorptivity and solar flux terms are as previously determined Solar flux is 
taken to be a fraction of that at the surface of the Earth (890 W/m 2 ). The maximum cross 
sectional area of the probe seen by the earth is DL. The percent of solar power reflected 
by the Earth (albedo) was estimated to be 39%. 

Typical values of the reflected power of the sun from the Earth on the probe are of 
the order 2. 84x10^ W at sea level and are approximately constant with altitude. 

F. PLANCK RADIATION 

1. Planck Radiation from the Atmosphere 

Planck radiation consists of three terms, that being absorption from the surround- 
ing atmosphere, that being absorbed from the Earth, and that being radiated away by the 
probe. It is represented by 

Qplanck Q A"^9 e _ QP- (34) 

Energy transfer in the form of Planck radiation from the atmosphere to the probe is 

9a — SA^P atmj (3 5) 

where: 

e= emissivity (0.065), 

A6= area of the probe presented to the atmosphere (m 2 ), 

P atm = thermal radiation emitted from the atmosphere (W/m 2 ). 

The emissivity term is taken as a constant and given the value of that for 300 K radiation. 
The area As is the total probe surface area. The radiation term P atm was estimated by inte- 
grating over the plot of radiance versus wavelength for sky radiation as seen from the 
Earth, given in [Ref. 29 :p. 3-75]. For a first order approximation the atmospheric flux on 
a sunlit summer day at 30°-40° latitude with an estimated probe field of view of 47t stera- 
dians was used. An approximate value for P atm of 250 W/m 2 resulted. 
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Typical values of radiation from the surrounding atmosphere on the probe are of 
the order 9.6x1 O' 7 W at sea level and are constant as a function of altitude. P atm is more 
correctly a function of the air temperature and therefore would vary with altitude, but 
these effects are small and were neglected in this work.. 

2. Planck Radiation from Earth 

Planck radiation from the Earth is given by 

q E =eA 5 Pc, (36) 

with: 

8= emissivity for 300K radiation, 

A 5 = area of probe presented to earth (m 2 ), 

Pe= radiation from Earth (W/m 2 ). 

Emissivity is as before and As is equal to the cross sectional area DL. The Earth-radiated 
flux was estimated using the same integration approximation technique as used for P atm 
above but now integrating over the upper portion of the radiance versus wavelength plot 
of the Earth as seen from space as given in [Ref. 29:p. 3-54], This integration yielded a 
first order approximation of the solar flux on a sunlit summer day at 30°-40° latitude with 
an estimated probe field of view of 27t steradians. An approximate value for P e of 615 
W/m 2 was used. 

Typical values of radiation on the probe are of the order 7.5xl0' 7 W at sea level 
and are constant as a function of altitude. 

3. Planck Radiation of the Probe Wire 

Planck radiation from the surface of the wire is 

q P =ecA 3 T 4 , (37) 

where: 

s= emissivity for 300K radiation, 
o= Stefan-Boltzman constant (5.67xl0‘ 8 W/m 2 K 4 ), 

A 3 = area of the wire probe radiating (m 2 ), 

T= atmospheric temperature (K). 
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The emissivity for platinum wire was 0.065 [Ref. 28:p. 10-282], The area term, A 3 , is the 
total surface area of the probe. The atmospheric temperature is used instead of the probe 
temperature T FO be to simplify the calculation and because AT«T probc as shown in [Ref. 
5:p 4]. 

Typical values of the Planck radiation emitted by the probe are of the order 
1.50x1 O ' 6 W at sea level and are shown as a function of altitude in Figure 3-20. 




Figure 3-20: Planck radiation from the 4.7 [im probe as a function 
of altitude with constant air velocity of 5 m/s. 



G. OHMIC HEATING OF THE PROBE WIRE 

Probe wire heating due to the current through the wire is 

( 38 ) 

where I is the current through the probe and R is the probe’s resistance. Brown in [Ref. 
14] measured this quantity to be 3.3x1 O ' 7 W at the small currents used in this type of 
probe, and we adopted his value for all of our modeling. 
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H. PROBE TEMPERATURE 



1. Steady State Free Stream Velocity 

Solving the equilibrium equation for Tprobe , using Eq. (15), the difference in probe 
temperature and atmospheric temperature AT versus altitude is shown in Figure 3-21. 




Probe Wire Differential Temperature (K) 



Figure 3-21 : Differential probe temperature versus altitude 
with constant free stream velocity of 0 and 5 m/s with all 
probe areas constant at their maximum value for a 4.7 pm plati- 
num coated tungsten probe. 

2. Normal Free Stream Velocity Fluctuations 
As the sonde rises, balloon pendulum motion and rotation occur. With the combined ver- 
tical horizontal and rotational velocity it is possible for the relative wind speed on the 
probe to drop to low values. The transient zero velocity condition is of interest in our 
model because forced convective heat transfer would approach zero. Near zero air veloc- 
ity, convection would be much smaller allowing the solar flux to dominate probe heating 
and increasing AT . In practice, though, free convection will create an induced effective 
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velocity over the probe wire, due to density differences near the wire so flow velocities 
will only approach zero. 

Figure 3-22 represents the effects of velocity fluctuations on the differential tem- 
perature as a function of altitude. 
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Probe Wire Differential Temperature (K) 

Figure 3-22. Simulated probe differential temperature for a 
normal velocity fluctuation of 2.5 m/s about a velocity of 5m/s 
versus altitude for the 4.7 pm probe. 



3. Variation of Solar Dlumination Area 

The effect on the differential probe temperature of varying the probe’s orientation 
with respect to the sun is shown in Figure 3-23 for sinusoidal change of area A 2 as a func- 
tion of differential temperature and with constant altitude. 
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Figure 3-23: Effect of varying the solar illumination area, A 2 , 
on the 4.7 pm probe differential temperature with constant al- 
titude and air velocity (5 m/s). 



4. Probe Orientation with Respect to the Earth 

By sinusoidal variation of the area terms A 4 and A 5 with periods of 5 and 20 sec- 
onds, the effect of relative Earth-probe orientations on the probe differential temperature 
for constant altitude and air velocity (5 m/s) may be calculated. Examples are shown in 
Figure 3-24. 
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Figure 3-24: Effect of sinusoidal variation of Earth illuminated probe 
areas A 4 and As, (i.e., reflected solar radiation and Planck radiation 
from Earth) on the 4.7 jam probe differential temperature with con- 
stant altitude and air velocity (5 m/s). 

I. SUMMARY 



Relative energy contributions to the heat flux of the platinum coated tungsten wire 
are ranked by order of magnitude, and hence contribution to probe temperature, in Table 
3-3 for an air velocity of 5 m/s at sea level. Probe convection and direct solar radiation 
terms dominate, with convection being strongly dependent on free stream air velocity. 
Solar radiation effects are strongly dependent upon probe orientation. 

Fluctuations in the free stream air velocity over the probe can produce a maximum 
change in probe differential temperature of 0.1 K at 10 km. 

Direct solar radiation loading alone accounts for a maximum probe temperature 
change of approximately 0. 1 1 1 K with constant free stream velocity, V=5 m/s at sea level. 
This change is a consequence only of a change in probe solar illumination area A 2 . 
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Variation of the probe cross sectional areas presented to Earth (A, and A 5 ) results 
in a maximum of 0.028 K variation of probe temperature with constant free stream veloc- 
ity of V=5 m/s at 30 km. 

Table 3-4 summarizes the maximum and minimum variation of the probe differen- 
tial temperature as a function of free stream air velocity, probe solar illumination area, and 
probe Earth illuminated areas. 



Heat Transfer Mechanism 


Power (pW) 


Convection 


-10.63 


Direct Solar Illumination 


7.28 


Earth Reflected Illumination 


2.84 


Probe Planck Radiation 


-1.50 


Atmospheric Planck Radiation 


0.96 


Earth Planck Radiation 


0.75 


Ohmic Heating 


0.30 



Table 3-3: Typical values for energy transfer mechanisms involved in the AFRL 4. 7 pm 
platinum coated tungsten probe wire at constant air velocity of 5 m/s at sea level with 
maximum wire areas illuminated. 



Parameter Varied 


Maximum 
Change in T 

(K) _____ 


Altitude 

(km) 


Minimum 
Change in T 
(K) 


Altitude 

(km) 


Free Stream Velocity 
Variation about 5 m/s 
(2.5a) 


0.1 


10 


0.036 


30 


Solar Illumination 
Area, A 2 


0.111 


30 


0.024 


0 


Earth Visible Areas, 
A, A 5 


0.028 


30 


0.011 


0 



Table 3-4: Maximum and Minimum change in values for the AFRL 4.7 pm probe differ- 
ential temperature for given parameter variation. 
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IV. EXPERIMENTAL MEASUREMENTS 



A. PROBE TEMPERATURE SENSITIVITY CHARACTERISTICS 

1. Purpose 

Verification of the actual probe response characteristics when subjected to solar 
radiation and varying wind speeds was done experimentally to validate probe modeling of 
the solar and convective heat transfer mechanisms. In this experiment we measured the 
differential temperature of the probe in response to one solar flux of radiation as a function 
of wind speed. The effect of Earth reflected sunlight at zero wind speed on the probe was 
also approximated experimentally. 

2. Experimental Procedure 

The platinum coated tungsten probe was mounted in a wind tunnel with air flow 
perpendicular to the wire (Figure 4-1). A Radiac ELH 300 W tungsten halogen lamp was 
focused uniformly on the probe through an optical lens, and the lamp power was adjusted 
to correspond to one solar flux after being chopped with a Laser Precision Corporation 
CTX-534 Radiation Chopper. Calibration of our unit solar flux source was done by direct 
comparison with the Sun at local solar noon on a clear day at latitude 40N, as measured 
by a United Detector Technologies model 255 silicon photodetector. A uniform spectral 
response measurement was made with a Scientec 362 Thermopile Power Energy meter to 
cross check the spectral response of the photodetector. To the accuracy of our measure- 
ments (about 1%), the halogen lamp flux and the solar spectral flux produced the same 
probe heating with both the silicon and broadband detectors. The signal was passed 
through a Stanford Research Systems SR560 Low-Noise Preamplifier and voltage read- 
ings were then taken using a Hewlett-Packard 3 561 A Dynamic Signal Analyzer and Ni- 
colet Pro 30 digital oscilloscope for wind speeds from 0 to 6.0 m/s. Wind velocities were 
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measured using an Omega HHF 52 hot wire anemometer accurate to 0. 1 m/s. Three sepa- 
rate experiments were conducted on 24-25 April, 16 May, and 29 May 1997 using three 
different platinum coated tungsten probes. Probes #1 and #2 were from similar batches of 
platinum coated tungsten wire while Probe #3 was a much older (1989) probe from a 
separate production run of wire obtained from the Air Force Geophysics Directorate. 

The effect of Earth reflected sunlight was measured using a 41% reflective gray 
surface behind the probe. The reflected sunlight measurements were made at zero wind 
speed. 



3. Results 

Experimental results are tabulated in Tables 4-1 to 4-5. The rms voltages read 
from the HP-3561 A spectrum analyzer were converted to peak-to-peak voltages by multi- 
plying by a factor of 2-jl and dividing by the gain applied at the SR560 preamplifier, 
since the wave form of the chopped light was nearly sinusoidal. For the second experi- 
ment (16 May 97) this procedure was modified slightly to include a correction factor, of 
2.22, to account for an observed non-square, chopped wave form from the halogen lamp 
source. The probe differential temperature was obtained by solving for AT in the equation 
for probe resistance [Ref. 13:p. 11] 

R w =R 0 {\ + aAT), (39) 

where Ro is the cold reference probe resistance. A different probe was used in each ex- 
periment so this value was measured as 19.94 Q (24-25 April 97), 18.81 Q. (17 May 97), 
and 25.68 Q (29 May 97). The constant alpha was assumed to be a=0.0042 for platinum 
coated tungsten wire [Ref. 30], 

An upward drift in measured voltages was noted after continual illumination of 
about 60 s or more. We suspect that this was most probably due to probe support heating 
effects as shown in Ref. 5:p. 16. Severe direct thermal contamination was noted with the 
halogen lamp closer than about 0.5 m, obviating the use of a lens to separate the hot lamp 
from the wind tunnel. 
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The observed temperature change depended on the chopping frequency and the 
probe response time. The probe response time varied between 0.6 to 2 ms and acted like 
a one pole low-pass filter at these frequencies. Collecting measurements at low frequen- 
cies, on the order of 5 Hz, was not always possible because of interference from ambient 
turbulence. In these cases higher frequencies around 25 Hz were used after including a 
correction factor. 

Probe #3 showed a lower probe differential temperature (Table 4-4) than Probes 
#1 and #2. Possible explanations for this discrepancy are investigated in Chapter V. 



Air Velocity 
(m/$) 


Measured Volt- 
age hp 3561 A 
(mV- rms) 


Peak-to-Peak 
Voltage - 
(mV) 


Actual Volt- 
age (nV) 


AR (mQ) 


Differential Tem- 
perature (K) 
[Corrected] 


0 


129.0 


364.9 


7.30 


7.97 


0.095 [0.1201 


1.0 


114.3 


323.3 


6.47 


7.06 


0.085 [0.091] 


2.0 


96.88 


274.0 


5.48 


5.98 


0.071 


3.0 


82.51 


233.4 


4.70 


5.13 


0.061 


4.0 


79.70 


225.4 


4.50 


4.91 


0.059 


5.0 


74.50 


210.7 


4.20 


4.58 


0.055 


6.0 


67.00 


189.5 


3.79 


4.14 


0.049 



Table 4-1: Experimental probe differential temperature measurements of Probe #1 (new) 
for one solar flux from a halogen lamp using a chopping frequency of 50 Hz and cor- 
rected temperatures for 0 and 1 m/s air velocity measurements (19 April 1997). 



Air Velocity 
(m/s) 


Measured Volt- 
age hp 3561A 
(mV- rms) 


Peak-to-Peak 
Voltage - 
(mV) 


Actual Volt- 
age (nV) 


AR (mQ) 


Differential Tem- 
perature (K) 
[Corrected! 


0 


128.1 


362.3 


7.3 


7.9 


0.095 [0.120] 


1.0 


114.2 


323.0 


6.5 


7.1 


0.084 [0.090] 


2.0 


98.6 


278.9 


5.6 


6.1 


0.073 


3.0 


88.6 


250.5 


5.0 


5.5 


0.065 


4.0 


81.5 


230.5 


4.6 


5.0 


0.060 


5.0 


76.8 


217.3 


4.4 


4.8 


0.057 


5.5 


73.8 


208.7 


4.2 


4.6 


0.054 



Table 4-2: -Experimental probe differential temperature measurements of Probe #1 (new) 
for one solar flux from a halogen lamp using a chopping frequency of 25 Hz and cor- 
rected temperatures for 0 and 1 m/s air velocity measurements (24 April 1997). 
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Air Velocity 
(m/s) 


Measured Voltage Nicolet 
Pro 30 (mV- rms) 


Actual Volt- 
age(jiV) 


AR (mQ) 


Differential Tempera- 
ture (K) 


0 


40.69 


9.19 


10.27 


0.130 


1.0 


25.85 


5.84 


6.52 


0.083 


2.0 


20.82 


4.70 


5.03 


0.067 


3.0 


18.56 


4.19 


4.68 


0.059 


4.0 


16.93 


3.82 


4.27 


0.054 



Table 4-3: Experimental probe differential temperature measurements of Probe #2 (new) 
for one solar flux from a halogen lamp using a chopping frequency of 5 Hz (16 May 
1997). 



Air Velocity 
(m/s) 


Measured Voltage Nicolet 
3091 (mV- rms) 


Actual Volt- 
age OiV) 


AR (mQ) 


Differential Tempera- 
ture (K) 


0 


37.2 


8.26 


9.19 


0.085 


1.0 


24.9 


5.53 


6.15 


0.057 


2.0 


20.3 


4.51 


5.02 


0.047 


3.0 


17.6 


3.91 


4.35 


0.040 


4.0 


16.0 


3.55 


3.95 


0.037 


5.0 


15.0 


3.34 


3.71 


0.034 


Table 4-4: ] 


Experimental probe dil 


fferential temperature measurements for Probe #3 



(1991-3.8 pm) one solar flux from a halogen lamp using a chopping frequency of 10 Hz 
(29 May 1997). 



Air Velocity (m/s) 


25 Hz, Measured 
Voltage (mV) 


5 Hz, Measured 
Voltage (mV) 


Correction Factor 


0 


32.14 


40.25 


1.252 


1 


24.79 


26.59 


1.073 



Table 4-5: Correction factor for low air velocity and higher frequency chopping rates. 



Adding a gray cardboard background for reflection of the halogen lamp provided 
the probe with a 27t steradian solid angle reflective background surface which roughly 
simulated that of the Earth. Table 4-6 shows the increased signal differences due to 
“Earth reflected solar radiation” from the absence of a reflective surface, gray cardboard, 
reference (white) surface, and a first surface mirror. 
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Background (reflectivity) 


Measured Voltage - HP 
3561A (mV- rms) 


Percent Change (%) 


None (0%) 


34.9 


0 


Gray (-41%) 


43.0 


19 


White (-100%) 


60.0 


42 


First Surface Mirror (-100%) 


60.8 


43 



Table 4-6: Experimental data for the effect of reflective backgrounds on differential probe 
temperature of Probe #2 (new) with one solar flux of direct radiation at zero air velocity, 
sea level, and chopping frequency of 5 Hz (29 May 1997). 



Sources of error for this experiment include the signal lost due to the high chop- 
ping rate in the earlier data. This is partly resolved with the correction factor in Table 4-5. 
Other minor sources of error involve probe support heating effects, inaccuracy in wind 
velocity measurements, and assumptions that the halogen lamp accurately reflects the 
spectral output of the sun. 

4. Summary 

At sea level, with no forced air flow over the probe, radiation heating countered by 
free convection and Planck radiation caused probe heating about 0.130 K above ambient 
temperature. This differential temperature decreased as forced air flow over the probe 
effected more efficient forced convection cooling to about 0.055 K at 5.0 m/s. 

Earth reflected solar radiation contributed an additional 19% increase in probe 
temperature at zero forced air velocity and would be expected to remain independent of 
air velocity and approximately independent of altitude. 
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B. MEASUREMENT OF PROBE WIRE RESPONSE TO DIRECT SOLAR 
RADIATION 



1. Purpose 

The probe response to direct solar illumination at sea level was measured under 
natural conditions to detect any additional error sources or procedural omissions in our 
laboratory measurements and simulations. 

2. Experimental Procedure 

Using a similar procedure to that described in Section IV-A, above, the platinum 
coated tungsten probe was subjected to direct solar illumination on the roof of the NPS 
physics building (Spanagel Hall) at approximately solar noon on a clear day. High back- 
ground turbulence forced the use of a spectrum analyzer with chopping frequencies of 25 
Hz (see Figure 4-2) but otherwise the experimental procedure followed that of Section 
IV-A. 



3. Results 

The results of this experiment are tabulated in Table 4-7. Multiple data points for 
each air speed where taken and averaged because of the difficult measurement conditions. 

Possible error sources include effects from heated air rising from the rooftop con- 
crete surface. This caused noticeable fluctuations as the fan pulled the hotter turbulent air 
over the probe. The wind speed on the roof varied from zero to almost 5 m/s. Although 
the experiment was sheltered from direct wind, some breeze was evident. Ambient at- 
mospheric turbulence made it difficult to collect probe solar illumination data accurately. 
A chopping frequency of 25 Hz facilitated in this but forced a correction to the results as 
discussed in the previous section. 
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4. Summary 

Probe temperature rises, as seen in Table 4-7, from direct solar illumination were 
about 30% lower than those for the laboratory measured halogen illumination, after ap- 
plying the correction because of the higher chopping frequency. The error was most likely 
introduced from ambient air flow fluctuations encountered on the rooftop. Considering 
these corrections and the unsteady air, agreement with experiment was good. 



Air Velocity 
(m/s) 


Measured 

Voltage 

(mV) 


Actual 

Voltage 

(UV) 


Rw (mfl) 


AT (K) 


Average AT 
(K) 

(Corrected] 


0 


25.45 


5.86 


6.54 


0.0789 


0.0782 




25.40 


5.84 


6.53 


0.0788 






24.30 


5.59 


6.25 


0.0753 


[0.0980] 




25.60 


5.89 


6.58 


0.0794 




1 


20.70 


4.76 


5.33 


0.0642 


0.0749 




25.80 


5.94 


6.64 


0.0805 






22.76 


5.24 


5.86 


0.0761 


[0.0804] 




25.40 


5.85 


6.54 


0.0788 




2 


17.21 


3.97 


4.43 


0.0540 


0.0534 




18.03 


4.15 


4,64 


0.0560 






16.22 


3.74 


4.17 


0.0503 




3 


15.97 


3.68 


4.11 


0.0496 


0.0462 




13.83 


3.19 


3.56 


0.0429 






14.18 


3.27 


3.65 


0.0440 






15.50 


3.57 


3.99 


0.0481 




4 


13.43 


3.10 


3.46 


0.0417 


0.0430 




13.84 


3.19 


3.57 


0.0430 






14.27 


3.27 


3.67 


0.0443 




5 


10.72 


2.47 


2.76 


0.0332 


0.0337 




11.63 


2.68 


3.00 


0.0361 






10.24 


2.36 


2.63 


0.0318 





Table 4-7: Data from rooftop experiment using a chopping frequency of 25 Hz (16 May 
97). 
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c. 



MEASUREMENT OF PROBE WIRE VIBRATION EFFECTS 



1. Purpose 

Historically, error sources in the balloon thermosonde data were thought to result 
partly from probe wire vibration effects. As air flow over the probe induces vibration in 
the probe fine wire filament, it would be subjected to periodically varying solar radiation. 
This vibration, along with the package pendulum motion and rotation, could induce a sig- 
nificant artifact in the data. 



2. Experimental Procedure 



With the probe mounted in the wind tunnel (Figure 4-3) the behavior of the thin 
wire filament was investigated with respect to air flow speed and orientational depend- 
ence on vibration. Illumination of the probe wire was by a Hughes 7.35mW helium-neon 
laser beam chopped at 90 Hz. The probe output signal was amplified by a Stanford Re- 
search Systems SR560 Low-Noise Preamplifier, passed through a Wavetek Brickwall 
bandpass filter (model 753A), and read on both a Hewlett-Packard 3561 A Dynamic Signal 
Analyzer and a Tektronix 2236 100 MHz oscilloscope. The bandpass filter was set at 80- 
100 Hz and the gain was optimized for the signal being read. The probe was also visually 
monitored using an American Optical microscope placed above the probe but out of the 
air flow. 

Turbulent flow over the wire was induced using a thin rod placed just to one side 
and upstream of the probe. Reynolds number for the flow over this rod is on the order of 
Re=1000, which is in the turbulent regime and corresponds to a vortex shedding frequency 
of 50 Hz from [Ref. 3 1 :p. 15] 




(40) 



with S the Strouhal number ( S=0.15 in this case), V the free steam velocity, and D the 
thin rod diameter. 
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3. 



Results 



Models for flow induced vibration of thin wires developed by Perry in [Ref 13 :p. 
99] and Blevins [Ref. 31:p. 15] predict vibration frequencies about two orders of magni- 
tude greater than what were observed visually. Because of the very low Reynolds num- 
bers involved with the platinum coated tungsten wire, the probe wire vibration is most 
likely due to turbulence in the free stream flow over the wire and not to resonance or 
vortex shedding effects from the wire itself. As previously shown, the Reynolds number 
throughout the realm of our model is less than 4 (Figure 3-6 and 3-7), which is under that 
for separation of flow and turbulence for a circular cylinder which occurs with Re>5 [Ref. 
21 :p. 52], Observed vibration of the wire was very dependent on the presence of turbu- 
lence in the air flow (necessitating the thin rod obstruction) and also showed a marked vi- 
bration onset near velocities of 2 m/s with oscillations continuing at nearly constant 
frequency and magnitude to the limit of our test, 5.8 m/s. 

With the probe oriented such that the wire was vibrating in and out of laser illumi- 
nation, the presence of filament vibration had a noticeable effect on the output signal char- 
acteristics as observed on the oscilloscope. Fluctuation about the mean sinusoidal 
(chopped) output would increase with visible turbulence but was sufficiently random and 
not of a magnitude that would allow spectral components or harmonics to be noted above 
the background noise on the Hewlett-Packard signal analyzer. 
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Figure 4-3: Laboratory bench setup for measurement of PL probe vibration effects. 






4. Summary 

Although the probe wire was noticeably vibrating in and out of the laser beam at a 
frequency of less than 50 Hz, near that of the thin rod vortex shedding frequency, no ap- 
preciable signal could be quantified to the accuracy of our experiment. From this we con- 
clude that probe solar illumination vibration effects are not a significant source of data 
contamination, although they can occur in the presence of velocity fluctuations in the air 
flowing over the probe. On the other hand, airflow over the probe will displace the probe 
wire, since it is not taught. Differential displacement of the probe wires, in a dual probe 
configuration, will alter the area exposed to the sun. The magnitude of this displacement 
depends on the drag force, which declines with altitude. This reduction with altitude can 
produce the pattern seen in Figure 2-5, where the solar induced contamination declines 
with altitude above 20 km. 

D. MEASUREMENT OF COEFFICIENT OF THERMAL CONDUCTIVITY 

1. Purpose 

Four analytical methods based on experimental data for the modeling of the coeffi- 
cient of thermal conductivity, h, have been presented in Chapter III. To provide a check 
of the analytical expressions we chose to independently measure h for our probe. Two 
experimental methods were used. Atmospheric pressure was varied from sea level to 1 
torr, with air velocities over the probe ranging from 0 to 5.8 m/s. Three different probes 
were measured with tungsten coated platinum wires from two different manufacturing 
batches. Probes #1 and #2 were purchased new for this work while Probe #3 was previ- 
ously used in 1991. 
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2 . 



Experimental Procedure 



The platinum coated tungsten probes’ response to a square wave input was meas- 
ured on both a Nicolet Pro 30 and Nicolet 3091 digital oscilloscope (Figure 4-4). Our 
first method used the heat rate of change equation given in Richards [Ref. 22:p. 87] 

^2 rj~> 

rp ~a = 2h ( T p rob ‘ - T )+ rk ~&r- 2c7£ »( r * ~ T o) » ( 41 ) 



where r is the probe radius, and T 0 is the effective black body temperature of the sur- 
roundings. Neglecting the last two terms in Eq. (42) (conduction and radiation terms) the 
estimated rise time of the probe response is: 



r = 



t£> 2 L 
V 4 . 



h(7£>L) 



pDc 

~4h~’ 



(42) 



where x is the rise time, c, is the heat capacity (133 J/kg-K), D is the probe wire diameter 
(m), and p is the density of tungsten (19300 kg/m 3 ) [Ref. 20, p. 543]. 

The second technique relied on measurement of the asymptotic temperature rise 
from Joule heating of the probe wire and the relations: 



q = -^- = I 2 R = hAAT, 

7?(T) = J?o(1-i-qA 7^, and 
V 2 R i 



R 



probe 



v x -v 2 • 



(43) 

(44) 

(45) 



Here V) is the voltage applied by the function generator with series resistance Ri, V 2 is the 
measured voltage output on the oscilloscope, and Rprobe is the probe resistance (Figure 4- 

5 ). 
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Figure 4 -4: Laboratory bench setup for measurement of PL probe coefficient of thermal conductivity, h 




Vi 



DC 

Voltage 




V2 



Figure 4-5: Circuitry for the measurement of probe resistance. 

Air velocity was varied from 0 to 5.8 m/s for both of the above methods and the 
minimum voltage possible was used to reduce overheat effects. Conduction effects 
through the probe wire supports were measured as described in Lomas [Ref. 19:p. 69] by 
placing the probe in a vacuum and measuring the heat transfer rate to the probe support 
needles by: 

Finally, one probe was subjected to a vacuum under a zero air flow condition and varying 
pressure as measured with a Vaisala PA 1 1 digital barometer. 

3. Results 

The three sets of data are shown in Tables 4-8 to 13. Values for h were calculated 
using the MATLAB program listed in Appendix D. 
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Sensitivity to probe overheating was found and is shown in Tables 4-13 and 4-14. 
This overheat condition was thus minimized as much as was possible when making meas- 
urements. The reason for this sensitivity is unknown and is still being investigated. 

Contamination of the probe wire through normal use resulted in an increase in the 
observed probe response time by about 20% over that of a thoroughly cleaned probe. 
This effect was thought to be negligible and was first discounted as an error source. Perry 
in [Ref. 13:p. 56] discusses its possible causes and prevention. 

Conduction of heat through the probe wire supports was measured to be 0.00 15W 
at an overheat of 53 K, zero wind speed, and a pressure of 0 mBar. This power flow de- 
creases as the square of the input voltage and is therefore very small for our voltages and 
corresponds to a 2. 8x1 O' 5 W/K temperature difference between the probe wire and the 
probe supports. 



Air Velocity 
(m/s) 


v L 

(mV) 


Vu 

(mV) 


V L +(l-l/e)AV 

(mV) 


1/e Rise Time 
(ms) 


h (W/m 2 -K) 


0 


448.08 


515.8 


490.89 


1.94 


1555 


1 


442.80 


491.2 


491.20 


1.14 


2646 


2 


443.60 


482.4 


467.96 


0.91 


3314 


3 


442.60 


476.6 


464.45 


0.75 


4022 


4 


442.80 


474.0 


462.52 


0.72 


4189 


5 


442.66 


471.2 


460.70 


0.65 


4640 


5.8 


442.77 


469.7 


459.82 


0.62 


4865 



Table 4-8: Experimental data for time constant measurement technique of Probe #1 (new) 
using Nicolet 3091 and varying air velocity(l May 97). 



Air Velocity 
(m/s) 


v L 

(mV) 


Vu 

(mV) 


V L +(l-l/e)AV 

(mV) 


1/e Rise Time 
(ms) 


h (W/m 2 -K) 


0 


362.5 


520.2 


462.2 


2.46 


1226 


1 


359.75 


461.75 


424.2 


1.19 


2535 


2 


360.5 


444.5 


413.6 


0.94 


3209 


3 


360.0 


432.0 


404.5 


0.87 


3467 


4 


360.0 


425.5 


401.4 


0.76 


3969 


5 


359.5 


421.5 


398.8 


0.73 


4132 



Table 4-9: Experimental Data for time constant measurement technique of Probe #2 (new) 
using Nicolet 3091 and varying air velocity (14 May 97). 
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Air Velocity 
(m/s) 


v L 

(mV) 


Vu 

(mV) 


V L +(l-l/e)AV 

(mV) 


1/e Rise Time 
(ms) 


h (W/m 2 -K) 


0 


450.4 


556.4 


517.4 


2.32 


1063 


1 


449.8 


510.6 




0.870 


2810 


2 


449.5 


499.2 




0.705 


3468 


3 


449.4 


492.2 




0.596 


4102 


4 


449.6 


489.8 




0.526 


4648 


5 


449.5 


485.8 




0.505 


4841 



Table 4-10: Experimenta 



Data for time constant measurement technique of Probe #3 



(1991) using Nicolet 3091 and varying air velocity (29 May 97). 



Air Velocity 
(m/s) 


V L (mV) 


Vu (mV) 


AT (K) 


h (W/m 2 -K) 


0 


584.0 


731.0 


96.7 


3184 


1 


573.8 


675.0 


• 


4582 


2 


573.5 


658.5 




5456 


3 


573.5 


647.5 




6267 


4 


573.5 


640.7 




6900 


5 


573.5 


637.8 




7210 


5.8 


573.5 


633.6 


37 


7712 



Table 4-11: Experimental Data for asymptotic temperature change technique of Probe #1 
(new) using Nicolet 3091 and varying air velocity (2 May 97). 



Air Velocity 
(m/s) 


Vl (mV) 


Vu (mV) 


AT (K) 


h (W/m 2 -K) 


0 


450.4 


556.4 


46.0 


3853 


1 


449.8 


510.6 


28.9 


5638 


2 


449.5 


499.2 


24.1 


6584 


3 


449.4 


492.2 


21.1 


7428 


4 


449.6 


489.8 


19.9 


7810 


5 


449.5 


485.8 


18.0 


8574 



Table 4-12: Experimental Data for asymptotic temperature change technique of Probe #3 
(1991) using Nicolet 3091 and varying air velocity (29 May 97). 
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Pressure (mb) @ 
Input voltage (V) 


V L (mV) 


Vu (mV) 


AT (K) 


h (W/m 2 -K) 


1019.5 (2V) 


571.5 


717.5 


97 


3144 


0 (200mV) 


112.0 


121.8 


27 




12 


111.7 


113.3 


4 


1924 


12 (600m V) 


168.1 


174.0 


10 


1128 


49.7 


168.0 


172.8 


8 


1379 


98.4 


168.0 


171.9 


6 


1691 


151 


168.5 


171.8 


5 


2005 


201 


168.1 


171.8 


5 


2005 


201 (800m V) 


224.4 


233.7 


12 


1201 



Table 4-13: Experimental Data for asymptotic temperature rise technique for Probe #1 
(new) using Nicolet Pro 30 and varying pressure and zero wind speed (5 May 97). 



Input Volt- 
age (V p-p) 


V L ( mV) 


Vu (mV) 


1/e Rise Time 
(ms) 


Overheat (K ) 


h (W/m 2 -K) 


1 


279.2 


286.8 


0.750 


10 


6475 


2 


259.8 


321.6 


0.725 


40 


8120 



Table 4-14: Sensitivity of coefficient of heat transfer to overheat condition (8 May 97). 



4. Summary 

The values for the coefficient of heat transfer, h, using the rise time method are 
consistent with those associated with the calculated ATs from Section IV- A. The values 
of h derived from the time constant method, Tables 4-8 and 4-10, are almost one half 
those from the asymptotic temperature rise technique. Tables 4-1 1 and 4-12. 

The values of h are very sensitive to the amount of overheat, and therefore cur- 
rents through the probe should be minimized to reduce error. Although measurements 
were made with four significant figures it was hard to achieve 20% accuracy in the values 
of h. 
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E. MEASUREMENT OF PROBE WIRE RESISTANCE 



1. Purpose 

The previous section measures resistance changes with probe current dynamically 
using digital oscilloscopes. Another method of determining h measured the probe resis- 
tance change with varying input current and changing wind velocity for static changes in 
the resistance. This version of the method in Section V-D relies on the accurate input of 
current to the probe and precise determination of current through the circuit using meas- 
ured voltage across the load resistor, Figure 4-5. 

2. Experimental Procedure 

As shown in Figure 4-6 the probe resistance was measured using a Fluke 8840A 
multimeter and a 5 l l 2 digit Keithley 197 Autoranging Microvolt DMM in response to a 
current input from the Hewlett Packard 623 5 A power supply. Different supply current 
levels were tried in order to minimize the overheat effects and change in load resistance 
with temperature. Relations for the calculation of AT and the coefficient of heat transfer 
have been developed in the previous section. Contact resistance in the probe was meas- 
ured and considered in the results. Accurate determination of the load resistance was 
done with three independent instruments using the four wire measurement technique [Ref. 
36], 
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Figure 4-6: Laboratory bench setup for measurement of PL probe resistance. 




3. 



Results 



The results using 

P = I 2 R = hAAT (47) 

are shown in Table 4-15. The contact resistance of the probe was measured as 0.488Q 
and has been subtracted from the tabular value of R 2 shown. 

Possible error sources include a thermocouple type response of the probe wire 
supports and some heating effects from the Hewlett-Packard 331210A function generator. 



Air Ve- 
locity 

(m/s) 


V! (Volts) 


Ri(Q) 


Ii (mA) 


V 2 (Volts) 


Ri(0) 


AT (K) 
with 

Ro=19.285 


h (W/m 2 K) 


0 


9.9847 


999.15 


9.9932 


0.203364 


20.3502 


13.1511 


2521.9 


1 


9.9889 


<c 


9.9974 


0.199086 


19.9138 


7.7633 


4183.9 


2 


9.9901 


cc 


9.9986 


0.197777 


19.7805 


6.1175 


5275.3 


3 


9.9909 


cc 


9.9994 


0.197103 


19.7115 


5.2656 


6108.3 


4 


9.9913 


cc 


9.9998 


0.196650 


19.6654 


4.6965 


6833.0 


5 


9.9915 


(( 


0.0100 


0.196344 


19.6344 


4.3137 


7428.0 



Table 4-15: Results o: 



probe resistance method for determination of the coefficient of heat 



transfer of Probe #2 (9 May 1997). 



4. Summary 

Values for this method, which is a refinement of that used for Table 4-11 and 4-12, 
are also higher than predicted from either our model or the rise time technique. 
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V. 



DISCUSSION 



A. MEASURMENT OF PROBE WIRE DIAMETER AND COMPOSITION 

1. Purpose 

All of the preceding experimental measurements and simulations rely on, among 
other parameters, the accurate knowledge of the platinum coated tungsten probe’s di- 
ameter. The manufacturer, Thermo Systems Incorporated (TSI), gives this as 3.8 pm, but 
data taken in the course of this research raised doubt as to accuracy of this number. 

Also questioned during the course of our research were the surface characteristics 
of the platinum coated tungsten wire as they relate to its reflectivity and heat transfer 
characteristics. 

2. Experimental Procedure 

A new sample purchased directly from TSI in 1997, and an older one obtained 
from the AFRL in 1991, were placed in the chamber of a Cambridge 20 KeV scanning 
electron microscope (SEM). The SEM had recently been calibrated and provided a scale 
accuracy of about 10 nm for our experiment. Measurements of the diameters were then 
taken at several positions along the wires, along with images of the regions measured. 
Images of the probe wires were also stored to disk for later analysis. A second sample of 
a Dantec pure platinum 5 pm wire was also imaged to provide further information about 
other vendor wire surface characteristics. 

A Kevex Energy Dispersive X-ray (EDX) spectrograph then took representative 
measurements of the elemental composition at locations of the newer wire. These were 
recorded both where the platinum coating was present and at locations where spalling, 
thermal cycling, and possibly other processes had removed large patches revealing the 
underlying bare tungsten. 
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3. 



Results 



Conservative measurements of the wire diameter at several locations, Figures 5-1 
and 5-2, for the new wire revealed an actual diameter of 4.7 fjm. Measurements of the 
older wire sample Figure 5- 3 showed its diameter to be 3.7 jim. 

The characteristics of the new wire revealed an average of up to 22% bare wire 
(absence of platinum coating) over most of the surface with several large patches of bare 
tungsten. Figure 5-4 and 5-5. Striations in the platinum coating, very similar to the bark 
on a redwood tree, created a very rough and non-uniform reflective and radiative surface. 
(Figures 5-l,2,4,5,and 6). 

In contrast to the rough surface seen in the new wire, Figures 5-1 and 5-2, the 
older wire’s surface is significantly smoother and more reflective, Figure 5-3. Although 
the image quality is not quite as good as those for the newer wire, the lack of deep stria- 
tions and dark patches is still discernible. The much smoother surface texture of a Dantec 
PI 1 pure platinum probe is shown in Figure 5-7. 

Figure 5-8 shows that the tungsten areas were relatively pure in composition, while 
Figure 5-9 shows the significant amounts of copper that were detected by the EDX spec- 
trum analyzer in the platinum areas. This probably indicates the presence of either con- 
tamination or impurities in the wire material, possibly introduced during the 
manufacturing process. 
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Figure 5-1: SEM image of the AFRL probe wire (new) showing a conservative 
diameter of 4.66p.m. 




Figure 5-2: SEM image of portion of new wire with 4.75pm diameter and sig- 
nificant flaking . 
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Figure 5-3: SEM image of portion of older (1991) wire with 3.73 pm di- 
ameter. 




Figure 5-4: SEM image of large bare tungsten region with striations (new). 
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Figure 5-6: SEM image of wire end. 
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Figure 5-7: SEM image of a Dantec PI 1 5 pm pure platinum wire surface. 
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Figure 5-8: X-ray emission spectral components of bare tungsten wire. 
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Figure 5-9: X-ray emission spectral components of platinum wire coating. 
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4. Summary 

The new sample of the TSI platinum coated tungsten probe wire diameter was 
measured to be 24% greater than quoted by its manufacturer, while the older wire ob- 
tained from AFRL was 3% smaller than stated. Surface characteristics of the new wire 
showed a darker surface than appropriate for pure platinum and thus will introduce an er- 
ror in the values for the emissivity and albedo assumed in the solar heating model calcula- 
tions. Other heat transfer related effects from the rough and flaked surface of the new 
wire batch cannot be determined within the scope of this thesis research. 

Overall, the older platinum coated tungsten wire showed very similar surface char- 
acteristics to the pure platinum wire, as expected. The newer wire was significantly differ- 
ent and showed marked surface roughness, bare patches, striations, and significantly larger 
diameter than specified. 

B. COMPARISON OF MEASURED AND MODELED COEFFICIENTS OF 
HEAT TRANSFER, h 

Figures 5-10 and 5-11 compare the coefficient of heat transfer as modeled in the 
MATLAB program of Chapter III, for Probes #1 and #3 and the two different measure- 
ment techniques, i.e., time constant and asymptotic temperature change. The large dis- 
crepancy using the asymptotic temperature change technique and the time constant 
procedure for both probe wires is not understood and is still under investigation. 
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Figure 5-10: Modeled and measured h, for the 3.8 pm Probe #3 for 
measured rise time and asymptotic temperature increase techniques. 




Figure 5-11: Modeled and measured h, for the 4.7 pm Probe #1 for meas- 
ured rise time and asymptotic temperature increase techniques. 
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C. COMPARISON OF MEASURED AND MODELED DIFFERENTIAL 
TEMPERATURES 



Figures 5-12 and 13 compare the differential temperature rise from solar and oh- 
mic heating effects, measured experimentally in Section IV-A, to temperatures computed 
using the MATLAB model of Chapter III. Measured Probe #3 data , in Figure 5-12, 
compares favorably to the models for emissivities of 0.435 (platinum coated tungsten) and 
0.444 (pure tungsten). The 4.7 jam Probe #2, Figure 5-13, does not correlate as well most 
likely due to the unknown heat transfer and absorption effects caused by the poor surface 
characteristics previously mentioned. 




Figure 5-12: Modeled (dashed and dotted)and measured (solid) data for the 3.8 pm plati- 
num coated tungsten wire Probe #3 differential temperature from solar and ohmic heating. 
Emissivities of 0.334 (pure Pt), 0.435 (Pt-W), and 0.444 (pure W) are shown for the 
modeled curves. 
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Figure 5-13: Modeled (dashed and dotted) and measured (solid) data for the 4.7 pm 
platinum coated tungsten Probe #2 differential temperature from solar and ohmic heating. 
Emissivities of 0.435 (Pt-W) and 0.444 (pure W) are shown for the modeled curves. 



D. EFFECT OF PROBE WIRE DIAMETER VARIATION 



The MATLAB simulation program developed for Chapter III provides a way to 
investigate the effect of different probe wire diameters on the differential temperature. 
Changes in the probe wire diameter with altitude on the differential temperature are re- 
flected in Figure 5-14. Simply reducing the wire diameter will reduce the probe’s differ- 
ential temperature. For example, a reduction in wire diameter from 5 to 2 pm would 
result in a AT reduction from 0.16 to 0.07 K or a 64% improvement, at 30 km. Further 
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reduction of probe diameter below 2 pm into the more pronounced non-linear region 
would have a larger effect. 




Figure 5-14: Effect of varying the AFRL probe wire diameter for 
the platinum coated tungsten wire on differential temperature at con- 
stant altitude and air velocity of 5 m/s. 



E. EFFECT OF PROBE MATERIAL AND DIAMETER CHANGES 



In light of the very poor surface characteristics discovered for the TSI wire, other 
probe wire materials were considered. The dominant parameter affecting the probe wire’s 
solar absorption is the spectral emissivity. The spectral emissivity of several favorable re- 
placement materials where calculated using the following relation: 
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0 



(48) 



£ \ 



V 




0 



where a(v) is the absorptivity as a function of frequency and 




(49) 



which is the Planck blackbody spectrum [Ref. 32:p. 17], Curves for the spectral emissivity 
versus energy from values tabulated in the CRC Handbook are compared with the corre- 
sponding black body spectrum of the sun, a 6000 K body, in Figure 5-15. The black body 
spectrum was used as the weighting factor for Figure 5-16, which compares the integrand 
of Eq. (48) for the four metals. Completing the above integrals, over the range corre- 
sponding to the transparent region of the atmosphere, about 0 to 195 nm, we then arrive 
at Table 5-1 which lists the spectral emissivity for the five wire materials investigated. 

Revisiting the previous section and using the spectral emissivities of the new mate- 
rials we can plot the change in differential probe temperature with varying wire diameter at 
a constant altitude and air velocity of 5 m/s as shown in Figure 5-17. 
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Figure 5-15: Emissivity versus energy for some probe materials of interest 
with 6000 K black body for comparison. 




Figure 5-16: Weighted emissivity curve for different probe wire materials. 
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Figure 5-17: Effect of change in spectral emissivity on differential 
temperature and wire diameter at 30 km and 5 m/s. 



Wire Material 


Tungsten (W) 


Platinum (Pt) 


Silver (Ag) 


Aluminum (Al) 


Spectral Emis- 
sivity 


0.444 


0.334 


0.116 


0.073 



Table 5-1: 



Average spectral emissivities of some se 



ected metals for the solar spectrum. 



In Figure 5-17 it is evident that a change of wire material from the present plati- 
num coated tungsten to silver would improve wire solar absorption performance by 67% 
with both wire diameters at 5 pm and at an altitude of 30 km. Reduction of the wire di- 
ameter to 2 pm and change in material to silver would reduce the solar heating by 85% at 
30 km. Since Or 2 is proportional to AT 2 a factor of 6.7 reduction in the solar induced 
temperature would reduce Cj 2 error by a factor of 45. 
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VI. CONCLUSIONS 



In this thesis, the measured values for the solar heating of the 4.7 jam platinum 
coated tungsten probe varied from 0.055 K at 5 m/s to 0.12 K at 0 m/s and 0.034 K at 5 
m/s to 0.085 K at 0 m/s for the 3.8 pm probe, at sea level. These values were 12% to 
34% larger than predicted by the detailed model for the 4.7 pm probe. Much better 
agreement occurred for the 3.8 pm platinum coated tungsten probe. Measurements and 
model results agreed to within 1-7%. Direct solar heating of the probe is significant. 
Temperature fluctuations over 0. 175 K, at 30 km, can be induced for both 3.8 and 4.7 pm 
probes through the periodic modulation of the rotation and pendulum motion of the sonde 
package both due to the direct illumination (heating) of the probe but also because of the 
changing air flow (cooling) over the probe (Table 3-4). Notably, the maximum change in 
the probe’s differential temperature occurs at an intermediate altitude of 1 1 km with de- 
creasing differential magnitudes occurring above and below this altitude. This effect is 
due to solar heating and is best shown in Figures 3-21 and 3-22. 

The uncertainty in measured versus modeled values for the differential temperature 
of the 4.7 pm platinum coated tungsten probe in the atmosphere is most likely due to the 
complicated surface characteristics of the TSI wire. Convection, radiation, and air flow in 
and around this rough surface introduces unknowns which could easily be of sufficient 
magnitude to account for the measured differences. 

The determination of the value of the coefficient of heat transfer, h, is important in 
characterization of the probe’s ability to transfer energy and was therefore modeled and 
measured with great care. Final values of h calculated through the model and measured 
experimentally agreed to within 20% at air velocities of 0 m/s and 10% at air velocities 
near 5 m/s. Probe wire surface characteristics most likely account for the consistently 
higher differential temperatures, of about 0.002 K, measured experimentally for the 4.7 
pm probe. 
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Susceptibility to the direct solar heating effects can be substantially reduced by 
minimizing the probe wire diameter. Additional improvement to solar contamination ef- 
fects can be realized by changing the probe material, notably to silver or aluminum. 
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VH. RECOMMENDATIONS 



The best solution to solar contamination of data would be to perform only night 
data collection. As this is not always an option in all research and tactical situations other 
options are necessary. The TSI platinum coated tungsten wire measured has been shown 
to have very poor characteristics, from a solar absorption standpoint, and thus would be a 
prime candidate for replacement. This may have been a result of manufacturing defects. 
Options include pure platinum wires which present a much more reflective and thermally 
smooth surface or another material such as silver or aluminum. The diameter of the probe 
wire is also a major contribution to solar contamination and should be reduced. Similarly, 
the thermal heating of the probe supports must be minimized using highly reflective mate- 
rials. 

Replacement altogether of the probe assembly has been proposed. Some replace- 
ments include an acoustic sounder sensor or a laser/optical device. If the old design probe 
must be used for balloon sounding, fiirther research could be performed into stabilization 
techniques for the sonde to minimize or remove the geometry dependence of the balloon - 
sun system. For the near term, reduction in the probe wire diameter and shift to a new 
wire material, such as silver or aluminum, is highly recommended. 
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APPENDIX A- AIRBORNE LASER PROGRAM SUMMARY 



Table A-l: Highlights 


of the Airborne Laser Program acquisition process [Ref. 6:p. 12]. 


September 1992 


Boeing and Lockheed receive separate contracts to assess how well existing large air- 
planes, the 707, 747, 767 and B-52 would perform while carrying a high-energy laser 
and beam control system. Both teams separately reach the same conclusion and rec- 
ommend a 747 platform to the Air Force as the best option for the mission. 


August 1993 


Team ABL Proposal submitted 


April 1994 (MSO) 


Team ABL and Rockwell selected as two competing teams for Concept Study Design 
($21 million, 33-month contract) 


March-November 1 995 


Boeing conducts series of “Flex Blue” tests to determine vibration environment for the 
Airborne Laser. 


April 1995 


TRW begins construction of BDL-2, (building block COIL). 


May-October 1995 


Boemg Wind Tunnel Testmg of Turret Designs. 


July 1995 


TRW TRACE facility produces first laser light to demonstrate 10 kW laser power goal 


September 1995 


Air Force Operational Concept Review (OCR). 


October 1995-present 


Lockheed Martin beam control tracking, pointing and adaptive optics compensation 
demonstrations 


January 1 996 


Draft RFP for PDRR phase released 


February 1996 


Lockheed Martin beam control system brassboard operational 


March 5-8, 1996 


Team ABL Concept Design Review with Air Force; Final Package due within 30 days. 


March 1996 


Second Draft RFP for PDRR released. 


May 1, 1996 


U.S. Air Force releases Airborne Laser final PDRR RFP. 


June 1996 


TRW plans demonstration of 200 kW Laser modules. 


July 1996 


Team ABL submits Airborne Laser PDRR Proposal Response. 


October 1996 


Final technical program update to Air Force. 


November 12, 1996 (MSI) 


Air Force selects Team ABL as one contractor for Airborne Laser Dem/Val. $1.1 billion 
5 year Program Definition and Risk Reduction (PDRR)contract to Team ABL 
(Boemg/Lockheed Martin/TRW) (cl) 


1996-2000 


PDRR Phase I (build one PDRR aircraft). 


2000-2001/2 


PDRR Phase II (test and fly PDRR aircraft. 


2001/2-2004 (MSD) 


Engineering Manufacturing and Development (EMD). 


2004-2008 (MSIE) 


Production (IOC 3 aircraft 2005/6, FOC 7 aircraft 2007/8). 
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Figure A-2: Airborne Laser Program acquisition timeline. 
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APPENDIX B - PROBE SIMULATION PROGRAM 



%Daniel J. Richardson 
%13 May 97 

%Probe Temperature vs. Air Temperature 
%Using Free and Forced Regimes 
%and Low Re number Equations 



%Set Constants: 
aipha=0 .435; 
epsilon=0 . 065; 
gamma=0 . 39; 
sigma= 5.67e-8; 
Pe=615; 

Pa=250; 
k=0 . 3 ; 

D=4 . 7 e-6; 

L=4 . Oe-3 ; 
v=5; 

betal=l . 458e-6; 
S=110 . 4 ; 
z 0=8.1; 
p0=101300; 

M=0. 0289644; 

R= 8 . 3 1 ; 
Psurf=890; 

W= 2.5; 
c=1010; 
g=9. 82; 
beta2=3 . 66e3; 



%Solar Absorptivity of Probe 

%Emmisivity for 300K Radiation 

%Solar Energy Reflected by the Earth (Albedo) 

%Stef an-Boltzman Constant [W/m A 2 . k~4 ] 

%Long Wave Radiation from Earth [W/m~2] 

%Long Wave Radiation Emmited by the Atmosphere [W/m^] 
%Atmospheric Transmission Term 
%Diameter of Probe Wire [m] 

%Length of Probe Wire [m] 

%Balloon Ascent rate 

%Dynamic Viscosity constant [ kg/smk~l/2 ] 

%Sutherland 1 s Constant [Kj 
%Scale Height [ km] 

%Scale Pressure [Pa] 

%Mean Molecular Weight [kg/mol] _ j 

%Universal Gas constant [J/mole K] 

%Solar flux at the surface of the Earth [W/m A 2] 

%Wind Velocity Factor 
%Specific Heat [J/kg K] 

%Gravity 

%Coefficient oc thermal expansion [1/K] 



%Atmospheric Variables: 
Altitude=0 : 10 : 30e3; 
kt=Altitude; 
Ps=Altitude; 
A2=Altitude; 
t=Altitude; 
?r=Altitude; 
Re=Altitude; 



%Altitude in [m] 
%Define array size 
%Define array size 
%Define array size 
%Define array size 



for i=l : length (Altitude ) ; 
if (Altitude ( i ) <11000) 

T ( i ) = (Altitude (i)/1000-44.62)/ (-0.15493); %Temperature of Atm. [K] (0-llkm) 

kt ( i ) = ( (Altitude ( i ) /1000 ) -50 . 57 )/( -1 983 . 47 ) ; %Thermal conduct. [W/m-K] (0-12km) 



elseif (Altitude ( i ) >=11000 & Altitude ( i ) <20000 ) 

T ( i ) *217 ; % 

kt ( i ) =0 . 01 95; % 



(ll-20km) 

(12-20km) 



else 

T (i) = (Altitude ( i) /1000 + 1 42. 75)/ 0.75; 
kt (i) = (Altitude (i) /1 00 0+206. 5) / 11636 ; 

end 

end 



% " (20-30km) 

% " (20-30km) 



Ps= (Altitude . "0 . 584 ) +890; 
P=p0*exp ( (-Altitude/ 1000) /z0 ) ; 
mu=beta 1* ( T . ^ (3/2) ) ./ (T+S) ; 
rho= (M.*P) . / (T. *R) ; 
nu=mu . /rho; 



%Solar Flux [W/m~2] 

%Pressure [Pa] 

%Dynamic Viscosity 
% Density 

%Kinematic viscosity [m~2/s] 



%v=W*randn (size (Altitude) )+5; 



%Random Air velocity generation around probe 



deltaT=0 . 1 ; 

Gr= ( ( rho. ~2 ) *g*beta2*deltaT* ( D A 3 ) ) . /mu; 
Re=D*abs ( v) . /nu; 

Reef=(Re.~2+ (0. 9*Gr.~0.418) .~2) ."0.5; 
Pr= (mu. *c) . /kt; 



%Grashof number 
%Reynolds Number 
%Effective Reynolds 
%Prandtl Number 



number mix . Convection 
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%Check for Free Convection ONLY regime 
%p.!3l Bryan 



for j=l : length (Altitude ) 
if ( Re ( j ) <0 . 0132 ) 

h (j )=(0.32*kt (j ) ./D) ; 



else 

h ( j ) = ( 0 . 24+0. 56*Reef ( j ) . ~0 . 45 ) . * ( kt ( j ) . / D ) ; %Mixed convection 

end 

end 



Nu=D*h./kt; 



%Nusselt Number 



%Areas Seen by the 
Al=pi* D*L; 

A2= ( D*L ) ; 
t=0:0. 25:750; 
ps=randn ( 1 ) ; 
pe=randn ( 1 ) ; 
wp= ( 2*pi ) /20; 
wr= ( 2*pi ) /5; 



Probe : 

%Surface Area of Probe wire 
%Area of Probe seen by Sun 

%phase factor for sun area 
%phase factor for the Earth area 
%angular frequency of Pendulum motion 
%angular frequency of Rotation 



A2=abs ( 1 . 52e-8* ( sin (wp»t-rps ) +sin (wr*t+pe ) ) ) /2 ; %Sinusoidal Sun Area Variation 

A4=abs ( 1 . 52e-8* (sin ( wp*t+ps ) +sin (wr*t+pe ) ) ) /2; %Sinusoidal Earth Area Variation 



A3=A1 ; 

A4= ( D* L ) ; 
A5=A4 ; 
A6=A1 ; 



%Area of Probe Radiating 
%Area of Probe seen by Earth 

%Area of Probe seen by Earth reflected Sunlight 
%Area of Probe seen by Atmosphere 



%Energy to/ from the Probe: 

Solar= ( alpha*A2 . * Ps ) ; 

Planck= ( eps i Ion* sigma* A3* T. ^4 ) ; 
EarthS= ( alpha*A4 * Psurf*gamma ) ; 
EarthP= ( k*epsilon*A5* Pe ) ; 

AtmP= ( epsilon*A6* Pa ) ; 



%Energy from the Sun 

%Planck radiation term for Probe 

%Reflected Energy of Sun from Earth 

%Long Wave Planck Radiation from the Earth 

%Long wave Planck Radiation from Atmosphere 



Tprobe= ( Solar-Planck+EarthS+EarthP+AtmP+Al*h . *T ) . / ( h . *A1 ) ; %Final Probe Temp 
%deltaT= . 095; 

Convection=h . *A1 . * (Tprobe-T ) ; %Convective Heat Transfer from wire 



%semilogx (Tprobe-T, Altitude/ 1000 ) ; 
plot(h, Altitude/1000) 

%grid; 

%title (' Probe Delta Temperature vs. Air Temperature' ) ; 
ylabel ( ’Altitude (km)'); 
xlabel ( ' h ' ) ; 
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APPENDIX C - DATA ANALYSIS PROGRAM 



MATLAB Data Analysis program developed by Prof. Donald L. Walters 

% vizdata.m program to read the Sept 90 Anderson peak 
% 

% 21 Apr 1997 : DLW 
% set up constants 
fsample=l/l . 25; 
f0=l/ (2*pi) ; %f ilter 

nf ft=256; 
n2=nf ft/ 2 ; 
n=52 ; 



% setup up file parameters 

fid=fopen ( ' c:\Met data\Ap092716 . txt ' , * r ' ) ; 
col s=7 ; % # of columns in the data set 

col=5; % column # of temperature data 

% read the file of arbitrary size 
[a, count] =f scan f (fid, ’ % f\t ’ ) ; 
status=f close ( f id) ; 
rows=count/cols; 
a=reshape(a,cols, rows ) ; 
a=transpose ( a ) ; 

% extract the high speed temperature data column 
b=a ( : , col ) ; 

c=b ( rows- ( nf ft-1 ) : rows ) ; % select end of array 

c=detrend (c) ; 
a=if ft ( c ) ; 
e=conj (d ) . *d; 

ps=e ( 2 : (n2 + l ) ) ; %subset of psd without DC term 

% create frequency array 
f=l:n2; 

f=f * fsample/nf ft ; 
plot ( f ,ps ) ; 

% create RC high pass filter 
fl=(f/f0) ."2; 

f ltr=fl . / ( 1 + fl ) ; %One pole filter 

f ltr=f ltr . ~f ltr ; %Filter applied twice 

% spectral variance fitr * psd (with + and - frequencies) 
ps2=ps . *f ltr ' ; 
semilogy ( f , ps2 ) 
var=sqrt I 2 . *sum ( ps2 ) ) 

% time domain filter 

[ B, A] =butter ( 1 , 0 . 16/ . 4 , ' high ' ) ; 
b2=filter(B,A, b-b ( 1 ) } ; 

b3=f ilter ( B, A, b2 ) ; 
var2=std (b3 (100: rows ) 5 

% Compute Cn2 

z=a(:,2); % Altitude 

cn2=b3 .'*’b3.* (79e-6* (a ( : , 3) . / (a ( : , 4 ) . *2 ) ) ) . "2; 

% Lowpass smooth the data 
( D, C] =butter ( 2 , . 05 ) ; 

% cn2=f ilter ( D, C, cn2 ) ; 
semilogx ( cn2 , z, 1 r ’ ) 
axis ( ( 1 . 0e-20 1.0e-13 0 20000]) 



% create one pole RC 
% apply filter 
% apply filter again 
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APPENDIX D- COEFFICIENT OF HEAT TRANSFER PROGRAM 



%Daniel J. Richardson 
% 4 May 97 

%Coefficient of Heat Transfer Measurement 



%Define Constants 
Alpha=0 . 0042 ; 

L— 4 . l5e-3; 

D=4 . 7e - 6 ; 

Rl=50; 

R0=G .485; 

Vl=2 - 0; 



%Constant for Platinum Coated Tungsten Wire 
% Probe Wire Length 
% Probe Wire Diameter 
%Resistance across Source 
% Probe Lead Resistance 

%Input Voltage form Function Generator 



%Measured Voltages 

V2=0,2244; %Voltage Measured at Low Point of Rise 

V3=0 .2337 ; %Voltage Measured at Asymptotic High 



%Calculate Required Terms 
A=pi*D*L; 

R2=(R1*V2}/ (V1-V2J-R0; 

R3= ( R1 *V3 ) / ( V1-V3 ) -R0; 
dT= ( R3-R2 ) / ( R2*Alpha ) 



%Area Of Probe Wire 
%Probe lead Resistances 
% " 

%Di f f erential Probe Temperature 



h= (V3^2 ) / ( A*dT*R3 ) 
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